Origin of the Temperature Oscillation in Turbulent Thermal Convection 
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We report an experimental study ol the three-dimensional spatial structure of the low-frequency 
temperature oscillations in a cylindrical Rayleigh-Benard convection cell. Through simultaneous 
multipoint temperature measurements it is found that, contrary to the widely accepted scenario, 
thermal plumes are emitted neither periodically nor alternatively, but randomly and continuously, 
from the top and bottom plates. We further identify a new flow mode — the sloshing mode of the 
large scale circulation (LSC). This sloshing mode, together with the torsional mode of the LSC, are 
found to be the origin of the oscillation of the temperature field. 

PACS numbers: 47.27.-i, 44.25. +f, 05.65.+b 



Thermal convection is a phenomenon that occurs 
widely in nature and in many industrial processes [l|. A 
paradigm to study the generic convection phenomenon 
is the Rayleigh-Benard (RB) system A fascinating 
feature of turbulent RB convection is the emergence of 
a well-defined coherent oscillation in the presence of tur- 
bulent background. This robust oscillation has been ob- 
served in both the temperature Q and velocity fields H, 
and in convection systems with different fluids [E [Ell 
and different geometries 0, B [§]• Although much effort 
has been devoted to the study of this phenomenon in 
the past [E EE [HI , the nature of this oscillation remains 
unsettled. It has been proposed that this oscillation is 
due to the periodic emission of thermal plumes from the 
top and bottom boundary layers of the system, which 
are coupled by the large-scale circulation (LSC) (hJ. In 
this scenario, plume emission due to boundary layer in- 
stability in one plate is triggered by the arrival of thermal 
plumes from the other plate, which implies that plumes 
are emitted not only periodically but also alternatively 
between the top and bottom plates. Some later experi- 
mental studies appear to support this picture and peri- 
odic plume emission has since been attributed to be the 
source of the observed temperature and velocity oscilla- 
tions in the bulk fluid 0, El, 0, El . In addition to the 
apparent oscillations observed for temperature and veloc- 
ity in the vertical circulation plane of the LSC, horizontal 
oscillations of the velocity field have also been observed 
HUES EH and it has been suggested that periodic plume 
emission is not necessary for the horizontal oscillation of 
the bulk fluid [IE El- Recently, it is conjectured that the 
local temperature oscillations may be caused by the tor- 
sional motion of the LSC 0]. The conjecture, which could 
explain the oscillation near the top and bottom plate, is 
unable to explain that observed at the mid-height plane, 
since, by symmetry, the torsional oscillations cancel out 
at the mid-height plane. We also note that the experi- 
mental studies that appear to show evidence of periodic 
plume emissions are two-dimensional (2D) measurements 
[U, EE El, El ■ With the LSC's twisting oscillation near 
the top and bottom plates and its azimuthal meander- 
ing [IE EE El, three-dimensional (3D) measurements 
become essential to unlock the intricate flow dynamics. 
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FIG. 1: (Color online) (a) Sketch of the convection cell and 
thermistor configurations. Power spectra of temperatures 
measured by thermistors (b) embedded in the top plate at 
positions TP1 to TP8 (from top to bottom); (c) at mid- height 
positions Tl to T8 in the bulk fluid; and (d) embedded in the 
bottom plate at positions BP1 to BP8. For clarity, each data 
set is shifted up from its lower neighbor by a factor of 20. 



In this Letter, we present 3D measurements of the tem- 
perature oscillations in a cylindrical convection cell of 
aspect ratio unity. Local temperatures are measured si- 
multaneously by 24 probes placed inside the convecting 
fluid and embedded in the top and bottom plates. These 
measurements reveal the origin of the temperature oscil- 
lation in the system and shed light on the driving mech- 
anism of its flow dynamics. The experiment was con- 
ducted in a RB convection cell that has been described 
elsewhere [la ]. Briefly it is a vertical cylinder of height 
H = 19.0 cm and diameter D = 19.0 cm, with upper and 
lower copper plates and Plexiglas sidewalk Water was 
used as the working fluid and measurements were made 
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at the Rayleigh number Ra = 1.7 xlO 9 and 5.0 xlO 9 , 
with Prandtl number Pr = 5.3. As the two measure- 
ments give the same qualitative results, only results for 
Ra = 5.0 xlO 9 will be presented. Local temperatures 
in the fluid are measured by 8 thermistors of 300 /zm in 
diameter. As shown in Fig. [Ha), these thermistors are 
mounted on a star-shaped frame made of 1 mm diameter 
stainless steel tube. They have an equal azimuthal sep- 
aration and have a distance of 1 cm from the sidewall, 
where the main ascending and descending flows are pass- 
ing through [TEl ]. The frame is soldered perpendicularly 
to a 1 mm diameter stainless steel tube and can thus 
traverse vertically along the central axis of the cell. The 
thermistors are labelled as 1, 2, ... 8, which also represent 
their azimuthal positions. The temperatures in the top 
and bottom plates were measured by 16 thermistors of 
2.5 mm in diameter [shown as solid circles in Fig. [I] (a)]. 
Eight of them are embedded in the top plate at half ra- 
dius from the plate center and about 2 mm away from the 
fluid contact interface at positions TP1 to TP8; while the 
other eight are similarly embedded in the bottom plate at 
positions BP1 to BP8. To lock the azimuthal orientation 
of the LSC steadily, we tilted the convection cell with 2° 
at position 1. Temperatures of the 8 small thermistors 
and 16 large thermistors are measured simultaneously by 
a 65-digit multimeter at a data rate of 0.55/sec. 

Figures W[b) and (d) plot the frequency power spectra 
of the temperatures measured by the embedded thermis- 
tors in the top and bottom plates, respectively. Due to 
the finite thermal diffusivity of the plates, when a hot 
(cold) plume is emitted from the bottom (top) plate, it 
leaves a cold (hot) spot there for a finite time before the 
temperature recovers to the previous value by conduc- 
tion. Similarly, when a hot (cold) plume impinge the 
top (bottom) plate, it will heat (cool) fluid in the bound- 
ary layer and leave a thermal imprint there. It has been 
shown in previous experiments that thermistors embed- 
ded in the plates are able to capture signatures of these 
plume departures and arrivals [ill [20| | . The prominent 
peaks in the power spectra near /o(— 0.028 Hz) corre- 
spond to the same oscillation frequency (for the same 
Ra) measured in many previous studies both in the fluid 
and in the plates [3, HI, EH E2, EH ■ It is seen from the 
figures that strong oscillations are observed at positions 
TP2, TP3, TP7, and TP8, where the hot ascending flow 
arrives at the top plate, and at positions BP3, BP4, BP6, 
and BP7, where the cold descending flow arrives at the 
bottom plate. On the other hand, no significant oscilla- 
tion is observed at positions TP4, TP5, and TP6, where 
cold plumes are emitted, and at positions BP1, BP2, and 
BP8, where hot plumes are emitted. These results show 
clearly that there is no temperature oscillation at po- 
sitions where plumes are emitted. Furthermore, if the 
plumes are not emitted periodically, then the oscillations 
observed at positions where plumes arrive from the op- 
posite plate must originate in the bulk. 

An important feature of the models that attribute the 
origin of temperature oscillations to the periodic emis- 
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FIG. 2: (Color online) Temperature cross-correlation func- 
tions measured in-plate at where plumes depart or arrive. 
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FIG. 3: (Color online) Segments of time traces of the temper- 
ature measured at the mid-height plane by thermistors T\ to 
T 8 . For clarity, Ti, T 2 and T 8 are shifted up by 1°C. 



sion of thermal plumes from the boundary layer is that 
emission of the hot (cold) plumes is triggered by arrival 
of cold (hot) plumes and thus hot and cold plumes are 
emitted alternatively from the bottom and top plates 
with a time separation equals to half of the oscilla- 
tion period. To test this, we examine the phase re- 
lationships between the temperatures measured in the 
plates, which can be obtained from the cross-correlation 
function of the measured temperatures: Ct^tAt) — 
{{Ti(t + T) - (T^iTjit) - {Tj)))/aT^ T] , where a T% and 
ax, are the standard deviations of the two quantities re- 
spectively. Figure [5] shows the cross-correlation function 
between temperatures measured at the positions BP1 
and TP5, where hot and cold plumes are emitted respec- 
tively. If the hot and cold plumes are shed alternatively 
from the bottom and the top plates, the cross-correlation 
function Cbpi,tps(t) should have strong negative peaks 
at r = ±To/2, with tq/2 the approximate time for plumes 
to traverse the height of the cell (here To = l//o = 36 
sec) . As shown in the figure the correlation is very weak 
and there is no negative peaks at r = ±r /2, but a nega- 
tive peak at r = 0. A negative peak at t = means that 
the hot and cold plumes arc shed simultaneously. Also 
plotted in Fig. [2] is the cross-correlation function between 
temperatures measured at BP1 and BP5, the positions 
where in the bottom plate hot plumes depart and cold 
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ones arrive, respectively. If the arrival of the cold plumes 
at BP5 triggers the emission of the hot plumes at posi- 
tion BP1, the cross-correlation function will have a strong 
positive peak near r = 0, which is clearly not the case as 
shown in the figure. 

Figure [TJc) shows the power spectra of the individual 
temperatures T\ to Tg measured in the bulk fluid at the 
mid-height (H/2) plane by the 8 small thermistors. The 
prominent peaks in the power spectra have the same fre- 
quency /o as in Figs, [H[b) and (d). It is seen that no 
significant oscillations are observed at positions 1 and 5, 
which are within the circulation plane of the LSC. On the 
other hand, strong oscillations are present at positions 2, 
4, 6 and 8. To understand why oscillation is present at 
some positions while absent in others at the mid-height 
plane, we plot in Fig. [3] segments of time traces of T x to 
T§. It can be seen that at positions 1, 2, 8 positive spikes 
dominate, suggesting hot plumes are ascending at these 
positions. At positions 4, 5 and 6 negative spikes dom- 
inate, suggesting cold plumes are descending there. At 
positions 3 and 7, both hot plumes and cold plumes are 
present. These temperature signals are consistent with a 
coherent large scale circulatory flow with a band of about 
half a diameter wide. To examine the phase relation- 
ships between the 8 temperature signals, we study their 
cross-correlation functions. Figure BJa) shows the cross- 
correlations between the (T2,Tg) pair and between the 
(T 4 , Tg) pair. The peak of Ct 2 ,t$ (Ct 4 ,t 6 ) near r = r /2 
indicates that the hot (cold) fluids pass the positions 2 
and 8 (4 and 6 ) alternatively with a time delay of tq/2. 
We calculate Ct^-t^ for a pair of thermistors when one 
of them mainly senses the upward spikes and the other 
mainly downward spikes. Figure BJb) plots Ct 2 ,-t 4 and 
Ct 8 ,-t 6 and the peak near r = indicates that the hot 
fluids pass 2 and cold ones pass 4 simultaneously. Similar 
relationship exists between positions 6 and 8. Simulta- 
neous presence of hot bursts at position 2 (8) and cold 
bursts at position 4 (6) is another evidence that hot and 
cold plumes are not emitted alternatively. The phase 
relationships above could also be observed from the tem- 
perature time traces shown in Fig. [3l Together, these 
results suggest a horizontal oscillation of the bulk fluid 
along the direction perpendicular to the vertical plane 
containing positions 1 and 5. In this off-center or slosh- 
ing mode of the bulk fluid, hot ascending fluid oscillates 
between positions 2 and 8 and cold descending fluid os- 
cillates between positions 4 and 6, and this explains why 
the oscillation strength at these positions are approxi- 
mately the same [Fig. 1(c)]. It should be noted that 
signatures of this sloshing mode has been observed by 
Qiu et al. [21] in a previous local velocity measurements, 
who found that the strongest oscillation occurs in the 
direction perpendicular to the LSC plane. 

To study the sloshing oscillation of the bulk fluid quan- 
titatively, we determine the hottest and coldest azimuthal 
positions of the bulk fluid along the sidewall from the in- 
stantaneous azimuthal temperature profile measured by 
the 8 thermistors at the mid-height plane. These po- 
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FIG. 4: (Color online) Temperature cross-correlation func- 
tions measured in-fluid at mid-height plane. 



sitions are determined by making a quadratic fit around 
the hottest and coldest temperature readings respectively 
[22| . The line connecting these two positions is the cen- 
tral line of the LSC band. The orientation of this line 
is the orientation of the LSC's vertical circulating plane, 
which is found to be very similar to that obtained by 
fitting a cosine function to the instantaneous azimuthal 
temperature profile [22| • We define the distance between 
this line and the cell's central vertical axis as the off- 
center distance d oc , which exhibits a periodic oscillation 
with an amplitude about 1/3 of the cell diameter [Fig. 
[UJa)]. Figure [5jb) shows the power spectrum of d oc , 
which has a peak at /o (0.028 Hz). Also plotted in the fig- 
ure is the power spectrum of d oc when the cell is levelled. 
It shows that the sloshing mode exists in both tilted and 
levelled cases and that it is not a result of the LSC being 
locked in a fixed azimuthal orientation. Figure^c) draws 
schematically this sloshing mode of the bulk fluid at mid- 
height plane based on the measured properties of d oc (t). 
This picture confirms some of the above conclusions in- 
ferred from the properties of cross-correlation functions 
between the local temperature probes. For example, it 
is seen that due to this periodical sloshing motion of the 
LSC, the thermistors at positions 2 and 8 (positions 4 and 
6) will alternatively experience the hot ascending (cold 
descending) flow. For the positions 1 and 5, as the LSC 
has a width of roughly half cell diameter, they will al- 
ways experience the hot ascending and cold descending 
flow but the degree of 'hot' and 'cold' varies due to the 
off-center oscillation. This explains the very weak peak 
at /o in the power spectra of T\ and T5. The peak at 
2/o of these spectra is due to the fact that within one pe- 
riod of the off-center motion the central line of the LSC 
crosses the (1,5) plane twice. From Fig. 5 (c) it is seen 
that the positions 3 and 7 are farthest from the band of 
LSC and so they can barely sense the bulk off-center os- 
cillation. This is evidenced by the barely visible peaks at 
/o in the power spectra of T3 and Tj. 

Some of the features observed by the thermistors em- 
bedded in plates can also be understood now. The os- 
cillation with frequency fo at positions TP2, TP3, TP7, 
TP8, BP3, BP4, BP6 and BP7, and the oscillation with 
frequency 2/o at positions TP1 and BP5 are due to the 
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the sloshing plus twisting motion of LSC. On the other 
hand, no significant oscillation is seen at positions BPf , 
BP2, BP8, TP4, TP5 and TP6. But when we placed 
the 8 small thermistors very near (3 mm) the bottom 
plates, positions 2 and 8 have strong oscillations. This is 
because that the horizontal positions of the plumes are 
modulated by the flow field only after they leave the re- 
spective plates and are thus not sensed by the embedded 
probes. This is a further evidence that the temperature 
oscillation does not originate from the boundary layers. 
It may be noted that some recent models of the LSC 
dynammics [HI . [23j ] have ignored the thermal boundary 
layer and the plumes. Thus, these models could in prin- 
ciple be modified to explain the results reported here. 
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FIG. 5: (Color online) (a) Segment of time trace of off-center 
distance d oc - (b) Power spectra of d oc for both levelled and 
tilted cells, (c) A schematic drawing of the off-center motion 
of the LSC at mid-height plane. The shaded band represents 
an instantaneous position of the horizontal cross-section of the 
LSC with its center at a distance d oc away from its average 
position (long-dashed line), and the two short-dashed lines 
represent the average position of the LSC band. The black 
dots represent the actual positions of the 8 small thermistors. 



With the above picture, we can understand the rea- 
son why some earlier measurements made in a 2D plane 
appear to observe that thermal plumes are emitted pe- 
riodically. In the present experiment the orientation of 
the LSC is locked steadily in the (1,5) plane by tilting 
the cell with a 2° angle so that it has a very small range 
of azimuthal angular fluctuation (<f> r ms — 8.5°). In some 
of the previous studies that observed temperature oscil- 
lations within the plane of the LSC, the convection cell 
was titled by lass than 1° [U, H ED]. When we tilted 
our cell by 0.5° it is found that the LSC is able to ex- 
plore much broader azimuthal range (<j) rms — 22.6°) [24J. 
When this azimuthal meandering is combined with the 
sloshing motion the LSC band could leave and return to 
the (1,5) plane so that the alternative occurrence of the 
hot and cold bursts in T\ and X5 can be observed. In fact 
in this case temperature oscillations at all 8 positions can 
be observed [22j| . Therefore the alternate appearance of 
hot and cold bursts at positions of 1 and 5 when the 
azimuthal orientation of the LSC is not steadily locked 
can be understood as a result of the horizontal motion 
of hot ascending and cold descending fluids being modu- 
lated by the sloshing mode. To conclude, our simultane- 
ous multipoint temperature measurements show directly 
and convincingly that thermal plumes are emitted nei- 
ther periodically nor alternatively from the top and bot- 
tom plates and that temperature oscillations are caused 
by the sloshing mode and the torsional mode of the ve- 
locity field in the central and boundary layer regions of 
the system, respectively. 
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